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Abstract 

Predictions on fluctuations of hadron production properties in central heavy ion collisions are 
presented. They are based on the Statistical Model of the Early Stage and extend previously 
published results by considering the strongly intensive measures of fluctuations. In several of the 
considered cases a significant change in collision energy dependence of calculated quantities as a 
result of the phase transition is predicted. This provides an opportunity to observe new signals of 
the onset of deconfinement in heavy ion collisions experiments. 
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1. Relativistic nucleus-nucleus (A-l-A) collisions provide a unique opportunity to study 
experimentally different phases of strongly interacting matter and transitions between them, 
for the recent review see Ref. [1]. In particular, since the discovery of sub-hadronic particles, 
quarks and gluons, it was expected that at high temperature and/or pressure densely packed 
hadrons will ’’dissolve” into a new phase of quasi-free quarks and gluons, the quark-gluon- 
plasma (QGP) |2H3]. 

Years of experimental and theoretical studies of high energy A-l-A collisions led to the 
conclusion that the QGP exists in nature. This conclusion is based on a wealth of system¬ 
atic data on A-l-A collisions at very high energies from the GERN Large Hadron Gollider 
(LHG) (see, e.g.. Ref. [6]) and the BNL Relativistic Heavy Ion Gollider (RHIG) (see, e.g.. 
Ref. 0 ). and, very importantly the evidence of the transition between hadronic matter and 
QGP (the onset of deconfinement) at the GERN Super Proton Synchrotron (SPS) ener¬ 
gies Pl9]. 

2. The experimental search for the onset of deconhnement was motivated [lO] by predic¬ 
tions of the Statistical Model of the Early Stage (SMES) [H] of A-l-A collisions. According 
to the model the onset of deconhnement in central A-l-A collisions should lead to a rapid 
change of the energy dependence of several hadron production properties, all appearing 
in a common energy domain. The predicted signals in single hadron properties were ob¬ 
served [HI El [12]. They indicate that the onset of deconhnement (the beginning of the mixed 
phase) is located at ^Jsnn{OD) ^ 8 GeV and the softest point (the end of the mixed phase 
region) at ^/s^{SP) ~ 12 GeV [13], where yG/jN denotes collision energy per nucleon pair 
in the center of mass system. However, up to now no convincing signal in huctuations of 
event properties was reported [TT] . 

Fluctuations are signihcantly more difficult to study than single hadron properties. This 
may explain problems in locating the onset of deconhnement using huctuation measurements 
and it asks for experimental and theoretical developments in the study of huctuations in 
relativistic nucleus-nucleus collisions. 

3. This work extends the model predictions [T^llG] on huctuations of hadron production 
properties in relativistic A-l-A collisions related to the phase transition. The predictions have 
been based on the SMES model of A-l-A collisions m- The extension is needed because new 
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measures of fluctuations were introduced in the recent years, for reviews see Refs. [HE]. 
In addition, the paper introduces a more general formalism to model fluctuations in high 
energy collisions which may be helpful in future efforts. 

4. Based on the success of statistical and hydrodynamical models of particle production 
in high energy collisions the SMES assumes that the matter created at the early stage 
of collisions is in equilibrium. The created matter has zero conserved charges, thus, all 
chemical potentials are equal to zero, and the hreball energy E and volume V are assumed 
to vary from collision to collision according to the probability distribution function P{E, V). 
Consequently, the energy density, £ = E/V , also may change from collision to collision. This 
leads to changes of other properties of matter which, within the grand canonical ensemble 
(GCE) used here, can be calculated using its equation of state (EoS). In particular, according 
to the first and the second principles of thermodynamics, the entropy change 6S {6 denotes a 
deviations from average value) is related to energy and volume changes a.s T6S = 6E+p6V, 
which provides T6S = V6e + {p + e)6V , where p is the pressure. Using the identity 
TS = E + pV one finds 

5S _ 1 5e 5V 

~S ~ 1+p/e 7 ^ ^ 

In case of energy and volume being proportional and the energy density being constant, 
5e = 0, Eq. Q gives: SS/S = 6V/V = 5E/E. Thus, the relative fluctuations of entropy are 
equal to those of energy and volume, and they are insensitive to the EoS. In other cases the 
entropy fluctuations depend on the EoS and thus on a form of the created matter. 

5. Similarly to entropy, mean multiplicity of particles of a given type changes with E 
and V. These changes depend on the EoS and particle properties providing they lead to 
changes of e. In particular, mean multiplicity of light pions or light quarks and massless 
gluons is approximately proportional to entropy of the produced matter. When crossing 
the transition region the effective number of degrees of freedom increases and thus the 
entropy increases faster with increasing ^/sNN■ Consequently, entropy fluctuations caused 
by energy density fluctuations will also be modihed. Mean multiplicity of strange hadrons 
in the confined phase or strange quarks in the deconfined phase is rather sensitive to their 
masses and effective number of degrees of freedom. These quantities are changed rapidly 
when crossing the transition region by increasing energy density. This results in a change of 
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fluctuations of mean multiplicity of strange particles caused by energy density fluctuations. 

6. Within the SMES the fireball created in a high energy collision is treated as a micro¬ 
state which belongs to the ensemble of all possible micro-states defined by the EoS and the 
energy-volume probability density function P{E, V). In particular, particle multiplicity is a 
property of a single event which varies from event to event even for fixed values of E and V. 

In the GCE these multiplicity fluctuations follow the Poisson distribution. Four extensive 
quantities E, V and N, Ns are considered. The distribution of the hrst pair, P{E,V), is 
assumed, whereas the distribution of the second pair V{N, Ns) follows from fluctuations of 
ensemble properties and fluctuations of event properties. 

Specially, the scaled variance of particle multiplicity can be presented [16] as a sum of 
the scaled variance induced by the initial E and V fluctuations, Uq, and the scaled variance 
of particle multiplicity at fixed E and V values, which is equal to unity for the assumed 
Poisson distribution: 

uj[N] = 1 + i:uo[iV] , uj[Ns] = 1 -|- uJo[Ns] . (2) 

7. For the purposes of this work it is sufficient to characterize the distribution P{E, V) 
by its five parameters which include its first and second moments: 

/A VW¥} (s V) - (£) {V) .... 

{ e ) • (V-) • y{(A£)2>((AV^)2> ’ ' ^ 

the latter three, two scaled dispersions and the correlation coefficient, are dimensionless. 

8. The SMES formulation and parameter values (unless otherwise stated) used in the 
original papers miiisiiisi are adopted here. Thus the ideal gas EoS is used to model the 
confined phase, the bag model EoS is used for the QGP phase, and the first order phase 
transition between them is assumed. 

9. Since event-by-event volume fluctuations cannot be eliminated in experimental studies 
of A-l-A collisions, it is important to minimize their effect by defining suitable fluctuation 
measures. It was shown within the model of independent sources that one can construct 
fluctuation measures from the first and second moments of two extensive event quantities, 
A and B, which are independent of the source number distribution. The measures were 
referred to as strongly intensive quantities na. The first measure of this type was introduced 
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in Ref. [IB] and then the concept was generalized na and extended to third [12] and higher 
moments PD]. Here the predictions of the SMES model will be calculated for strongly 
intensive quantities which include the hrst and second moments of A and B. 

10. Two families of strongly intensive quantities can be constructed HZI: 


Oa L 

E|AB] = - 2({AB) - {A)(B)) 


( 4 ) 

( 5 ) 


where oj[X] = ((X^) — (X)^)/(X). The normalization factors Ca and Cs are required to be 
proportional to hrst moments of any extensive quantities. Note that includes the 

correlation term {AB) — {A){B) whereas A[A,B] does not. 

11. Two selections of the Ca and Cs normalization factors are used in the present paper. 
Firstly, the normalization factors equal to mean of the second argument are assumed: 


Ca = Cj: = {B) . ( 6 ) 

As pointed out in Ref. pa this normalization within the statistical model of the ideal 
Boltzmann gas in the GCE (IB-GCE) formulation and R ~ R leads to: 

A[A,R] = S[A,R] = a;*[A] , (7) 

where ca*[A] is the scaled variance of A for a hxed system volume. 

Secondly, the normalization 


Ga = (X) - (Ns), Gs = (X) + {Ns) ( 8 ) 


will be used for particle multiplicities [2T]. It provides A[X, X 5 ] = S[X, X 5 ] = 1 in the 
IB-GGE. 

12. Within the IB-GGE and provided A and B are uncorrelated in a hxed volume, one 
hnds that uj*[A\ and uj*[B] can be expressed via S[A,R] and A[A, R]. Let us introduce the 
quantities: 


G[A,R] ^ - 


E[A,R] + A[A,R] = u[A] 


{A B) - (A) (R) 

(B) 


( 9 ) 
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^ i [E[B,yl] + A[B,yl]] = u,lB] - . 

Here the normalization of S and A is given by Eq. (§. Then one hnds [171120] 

a;*[H] = Q[A,B] , u*[B] = Q[B,A] . 

Predictions for f2 quantities will be also calculated in this paper. 

13. Four extensive event quantities E, V and N, Ns dehne their six pairs: 

\E,Vl l/VVI, lNs,V], IN,E], [Ns,E], [Ns.N], 


(10) 


( 11 ) 


( 12 ) 


for which the strongly intensive fluctuation measures are calculated. Within the model the 
results for the [E, V] pair depend only on the assumed distribution P{E, V). The remaining 
pairs include at least one extensive event quantity which fluctuations are dependent on the 
EoS. 

14. If the energy density £ remains constant from event to event, the energy and volume 

fluctuations are correlated as E ~ E. Within the IB-GCE these fluctuations do not influence 
the strongly intensive measures A and S. The fluctuations of e within the SMES lead, 
however, to their dependence on {V) and /{V). Namely they are proportional to 

the term 1 + ((51/)^)/(H)^. For central Pb+Pb collisions one gets 1 + (((5E)^)/(E)^ = 1 and 
thus in the following calculations the term is set to one. 

15. Based on the SMES model m the following numerical values of the parameters and 

their collision energy dependence are assumed. Mean volume and mean energy density at 
^/sNN = 10 GeV are set to be (V) = 350 fm^ and (e) =3.2 GeV/fm^, respectively. Their 
dependence on collision energy is taken to be {V) ~ 1/yGj^ and {e) ~ ^JsNN [yJsNN — 
2 ruAr). Scaled dispersion of volume fluctuations is set to / (V) = 0, which is a good 

approximation for central Pb+Pb collisions. Three values of the {{5eY) / {e) parameter are 
used: 0.17 (solid line), 0.13 (dashed line), and 0 (dotted line). The value ■>/((fe)^)/ {e) = 0.17 
is considered to be the upper limit based on the UrQMD and HSD simulations [22]. The 
correlation between E and V is set to zero. 

16. From Eqs. ( |4|[5|[^ one gets n[E,V] = A[E,V] = S[E,1/] = u:*[E] = {{5Ef)/{E) 

shown in Fig. as a function of for the model parameter values given above. This 

result is evidently insensitive to the EoS. 
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Figure 1: ^[E,V] = /\[E,V] = Yi[E^V] = u)*[E] = {{6E)'^)/{E) as a function of ^/snn in central 
Pb+Pb collisions at the CERN SPS energy range. See text for the numerical values of the model 
parameters and their y^slF/v-dependence. The normalization factors Ca = Cy, = (E) are used. The 
solid, dashed, and dotted lines correspond to ((5e)^)/(e) equal to 0.17, 0.13, and 0, respectively. 
The vertical lines indicate the beginning (the onset of deconhnement) and end (the softest point) 
of the mixed phase region, ^Jsnn{OD) = 7.4 GeV and ^Jsnn{SP) = 10.8 GeV, respectively. 



Figure 2; The same as in Fig. but for E] = A[iV, E] = S[A, E] = u;*[A] (left) and 
E] = A[A 5 , E] = S[iV 5 , E] = (right). The normalization factors Ca = Cy = (E) are 

used. 

The energy density fluctuations modify fluctuations of event properties like particle mul¬ 
tiplicity and these modifications are dependent on the EoS and particle type. On the other 
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Figure 3: The same as in Figs. and but for A[N,E] and T,[N,E] {left), and A[A^ 5 ,i?] and 
S[iV 5 ,F^] (right). The normalization factors Ca = Cs = (E) are used. 

hand the EoS and particle properties change signihcantly when crossing the phase transition 
region. Thus, one expects that the collision energy dependence of properly selected fluctua¬ 
tion measures may signal the transition region. The fluctuation measures which are sensitive 
to the EoS are plotted in Figs.[2]j^as functions of collision energy in the range which includes 
the phase transition region. Note that the N and Ns fluctuations for \/~{{SVy^/(V) = 0 
and ^((fe)2)/(e:) = 0 are assumed to be Poissonian and uncorrelated. 

17. The collision energy dependence of P] = A[iV, P] = S[A, P] = a;*[iV] and 
f2[iV5, P] = A[Ns, V] = S[A 5 , P] = a;*[A 5 ] is shown in Fig. |^for central Pb-|-Pb collisions 
at the CERN SPS energies. With ^^/((^e)^ / {e) decreasing to zero the results approach unity 
as expected in the IB-GCE. A and S are equal to each other for pairs of quantities [fV, P] 
and [A 5 ,P]. It results from the used approximations: N and Ns are uncorrelated in the 
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Figure 4: The same as in Figs. T][3 but for A[A^ 5 , N] (left) and N] {right). The normalization 

factors Ca = (-^) — {^s) and Cs = {N) + (Ns) are used. 


fixed volume, {N) ~ V and {Ns) ~ V, and {{^V)"^)/(y) = 0. The overall increasing trend 
of fl[N, V] and E] with y/s^N seen in Fig. [^is due to non-zero values of / (e). 

The modihcations of this energy dependence are observed in the region of the phase tran¬ 
sition. However, measurements of system volume (or quantity which is proportional to 
volume) are likely to be experimentally difficult or even impossible. 

Figure [^presents the fluctuation measures A[iV, E] and S[iV, E] as well as A[Al 5 , E] and 
S[A 5 ,i?] as a function of collision energy. Similarly to the results presented in Fig. |^the 
overall trend caused by the assumed energy density fluctuations is modihed in the region of 
the phase transition. The most pronounced modihcation is observed for A[A 5 ,F 1 ]. 

Finally, the collision energy dependence of A[A 5 , A] and S[A 5 ', A] is shown in Fig. 
The normalization (|^ is implied in this case. The most pronounced modihcation of the 
overall trend in the phase transition region is observed for A[A 5 , A], and a weak one for 
S[A5 ,A]. 

18. In summary, predictions on the collision energy dependence of huctuations of hadron 
production properties in central heavy ion collisions in the range of the phase transition are 
presented. They are based on the Statistical Model of the Early Stage and extend previously 
published results by predictions for the strongly intensive quantities. They are calculated 
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for six pairs of event quantities: 


[E,Vi [N,V], [Ns,V], [N,E], [Ns,E], [iV 5 ,iV], (13) 

where E and V stand for the system energy and volume, whereas N and Ns for multiplicities 
of all and strange particles. In several considered cases the collision energy dependence is 
signihcantly modihed in the phase transition region. The most pronounced change is seen 
for A[iV 5 ,iV]. This opens a possibility to observe signals of the onset of deconhnement in 
data on fluctuations of hadron production properties. 
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